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Summary 

An automated transonic design code has been 
developed which modifies the geometry of a wing 
from an initial shape in order to obtain a specified 
pressure distribution. This automated design code 
utilizes an existing design algorithm which relates 
changes in surface pressures to changes in the airfoil 
geometry. The design code iterates between this 
design algorithm and a three-dimensional transonic- 
small-disturbance analysis code. In addition to the 
wing, the analysis portion of the code can model 
a fuselage, pods, pylons, and a winglet, and this 
modeling allows the design to take place in a realistic 
flow environment. A two-dimensional option is also 
available for airfoil analysis and design. 

Several two- and three-dimensional test cases are 
described to illustrate the capabilities of the auto- 
mated design code. These test cases include repeats 
of two previous design efforts which did not use auto- 
mated computational design. The design test cases 
were generally successful in that the target pressure 
distributions were achieved in a reasonable number 
of design cycles. 

Introduction 

Recently there has been considerable progress in 
the development and application of transonic aerody- 
namic analysis codes. A number of codes have been 
developed that are capable of computing transonic 
flow around complex wing-body configurations. One 
particular analysis code, described in reference 1, 
solves the transonic-small-disturbance flow equation 
around a wing-body configuration. This code was 
later extended to include pods, pylons, winglets, and 
tails, and thus allows analysis of almost complete air- 
craft configurations (ref. 2). 

One important application of such analysis codes 
has been the design of airfoils and wings. Some 
examples of design goals include achieving a max- 
imum lift-drag ratio, a favorable pressure gradient 
to promote natural laminar flow, or a reduction in 
the strength of a shock. Depending on the situation, 
the wing (or airfoil) may be designed from an arbi- 
trary initial shape, or a known shape that exhibits 
desired characteristics may be refined to achieve the 
design goals. Typically, a starting shape is analyzed 
and the results are evaluated to determine what con- 
tour changes are needed to satisfy the design require- 
ments. The modified geometry is then analyzed to 
determine the effectiveness of the modifications. Al- 
though this is a cut-and-try approach, it has proven 
to be a successful design method (refs. 3 and 4). 
Several design methods have been developed which 
automate this approach for two-dimensional airfoils 


and three-dimensional wing sections (refs. 5 to 21). 
These algorithms have been formulated to eliminate 
shocks, to design to a specified pressure distribution, 
or to modify the initial geometry to minimize a spec- 
ified aerodynamic quantity. An immediate benefit 
of these automated design procedures is a reduction 
in man-hours necessary to complete the design. The 
feasibility of using such computational design meth- 
ods has increased as the accuracy and reliability of 
analysis codes have increased. 

In this study a transonic wing-airfoil design code 
has been developed which designs to a specified 
pressure distribution. A transonic-small-disturbance 
code (refs. 1 and 2) provides the necessary aerody- 
namic analysis of the initial and design geometries 
and allows the design to occur in the presence of a 
fuselage, pods, pylons, and winglets if desired. Be- 
cause of the geometric capabilities of the analysis 
code, more geometrically complex wing design prob- 
lems can be addressed than was previously possible. 

The remainder of this paper is organized in the 
following manner. A description of the analysis code 
is presented first, followed by a description of the 
design algorithm. Results from various two- and 
three-dimensional design test cases are then shown 
to demonstrate the capability of the design code. 
Finally, some comments and recommendations on the 
use of the design code are included. 
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Subscripts: 
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wing span 
pressure coefficient 
local chord 

function defining wing surface 
Mach number 
surface velocity 
iteration count 

physical Cartesian coordinates 
angle of attack 
design coefficients 
specific heat ratio 
incremental change 
2 y/b 

velocity potential 
subsonic relaxation factor 

indices indicating location of a 
point in the z-, 2/-, and ^-directions 


s value at surface 

x , y , z partial derivatives with respect to 

Cartesian coordinates 

oo free stream 

Description of the Aerodynamic Analysis 
Code 

The analysis code is described first since many 
features of the design method are dependent on the 
structure of the analysis code. The analysis code 
is the wing, body, pod, pylon, and winglet code of 
Boppe (refs. 1 and 2), which has been used to analyze 
a wide variety of configurations. A brief description 
of the analysis code is presented below with com- 
ments focused on a wing-alone configuration. Vis- 
cous effects, which may be included in the analysis 
calculations, are not discussed below. All analysis 
and design computations performed in this study are 
inviscid computations. Detailed descriptions of the 
code are given in references 1 and 2. 

The flow equation solved by the code is a non- 
conservative form of the transonic- small- disturbance 
approximation to the potential-flow equation: 


1 - Ml - 0 + l)Ml4>x ~ 


<f>XX 
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The classic small-disturbance equation is extended 
by retaining the term <t>\<t> xx to better approximate 
the velocity where the flow equation changes from 
elliptic to hyperbolic type. The terms <\> y <\> X y and 
<t> x <t>yy have also been retained to improve the res- 
olution of swept shocks. The solution is obtained 
by successive line overrelaxation (SLOR) on an em- 
bedded Cartesian grid system through the use of 
finite-difference techniques. A global, or crude, grid 
encompasses the entire flow field and extends to in- 
finity in all directions except in the y-direction, where 
symmetry about the x-z plane at y = 0 is assumed. 
Fine grids may be embedded in the global grid in 
regions of interest, such as around the wing, body, 
pod, pylon, and/or winglet (fig. 1). After an initial 
solution has been obtained on the crude grid (typi- 
cally after 80 to 100 iterations), a series of iterations 
are performed that alternate between the crude grid 
and the fine grid. The number of crude-fine itera- 
tions needed to obtain a converged solution varies 
between 200 and 600 depending upon the configura- 
tion and flow conditions. 


The crude grid in the x-direction is made up 
of 51 grid lines. Thirty-eight evenly spaced grid lines 
traverse the wing from the most forward portion to 
the most aft portion. Two tangent functions are used 
to stretch the remaining grid lines to upstream and 
downstream infinities. A portion of the crude grid is 
shown in figure 2. A tangent function is also used to 
generate the grid lines from the wing plane to plus 
and minus infinity in the ^-direction. In the spanwise 
direction an exponential function generates the grid 
lines such that 18 lie on the wing, with the wing tip 
located midway between the 18th and 19th grid lines. 
It is along these 18 grid lines that the wing fine grids 
are positioned (fig. 3). The boundary conditions for 
the symmetry plane are 

&-«} < 2) 

and for downstream infinity 

<t>yy + <t>zz = 0 (3) 

is enforced. The potentials for the remaining crude- 
grid boundaries are set to zero. 

The wing fine grids are evenly spaced in the 
chordwise and ^-directions. Grid points are located 
at every 1-percent-chord point beginning at 0.5 per- 
cent of the chord. This location of the grid points 
results in the leading and trailing edges being po- 
sitioned between grid points. The extent of these 
fine grids in the vertical direction can be determined 
by the user but is typically 30 percent of the aver- 
age wing chord above the wing and 10 percent below 
the wing. The wing fine grids extend 20 percent of 
the local chord in front of the wing and 10 percent 
behind it in the streamwise direction. The poten- 
tials along this outer boundary are determined by 
linear interpolation from the previous crude-grid so- 
lution and are held fixed. Neumann boundary condi- 
tions are applied at the wing plane to form an inner 
boundary. Once the wing fine-grid potentials are re- 
laxed, linear interpolation is used to compute crude- 
grid boundary points from fine-grid boundary points 
at the wing plane. These potentials then become 
Dirichlet boundary conditions for the crude-grid in- 
ner boundary. 

The wing is modeled along a constant -2 plane 
with double-valued potentials corresponding to the 
upper and lower surfaces of each wing section. In- 
formation concerning the shape of each wing section 
is transferred through the inviscid wing boundary 
conditions 


<t>z = fx~ a (4) 

where f{x,y) defines the wing surface. 


2 



Input to this analysis code includes the flow con- 
ditions (Mach number and angle of attack), wing 
planform definition, and wing section definitions. 
Spanwise and chordwise interpolations are performed 
on the input geometry to determine the 18 computa- 
tional wing sections. 

Description of Design Method 

The design method modifies a specified airfoil or 
wing geometry until a target pressure distribution is 
obtained. The design algorithm used in this method 
was developed by McFadden. (See ref. 13.) It de- 
termines the magnitude and direction of the changes 
to the wing (or airfoil) sections. The premise behind 
this algorithm is that at a given point the difference 
between the square of the target surface velocities 
and the square of the calculated surface velocities 
is proportional to a change in the surface ordinates. 
This equation can be expressed as 


2 2 __ R 

^target calculated “ ^0 


d{z s /c ) 
dt 


(5) 


terms are used regardless of the local Mach number. 
The modified method appears as 
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Using finite-difference techniques and allowing At to 
be unity (i.e., one design step), we may write the 
design algorithm for the j th point as 


\ f , ^ , * A (V c )j - A (zs/c)j-i 
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(8) 

Through combination of the A (z 3 /c) terms, the al- 
gorithm becomes 


where t is the iteration count and /?o is the propor- 
tionality constant. All changes made to the geometry 
are normal to the wing plane. (See fig. 4.) 

Two additional terms are added to this equation 
to obtain smooth design convergence. The impetus 
for selecting these two particular terms comes from 
the work of Davis (ref. 22). Davis observed that rela- 
tionships between the surface pressure coefficient and 
the first or second derivative of the surface (depend- 
ing upon the free-stream Mach number) can be de- 
rived from small-disturbance solutions to the wavy- 
wall problem (ref. 23). The resulting relationships 
are 

__ 2 dzs/dx 

P " (Ml -1)1/2 

for supersonic flow, and (6) 

_ £ d^Zg/dx^ 

P ~ 7T(1 — M2,) 1 / 2 

for subsonic flow, where t is the wavelength of the 
undulations. Based on these relationships for the 
pressure coefficient, Davis developed a predictor- 
corrector method which designs to a target pressure 
distribution. Surface corrections are determined by 
one of two design algorithms which utilize the wavy- 
wall relationships for the pressure coefficient. The 
choice of design algorithm depends on whether the 
local Mach number is less than or greater than unity. 
The two additional terms added to the McFadden al- 
gorithm (ref. 13) are the nondimensional time deriva- 
tives of the airfoil slope and curvature, and both 


A(zs/c)j +1 (^2 ) + A (W c )j (So 

+ A( z s/c)j-l ( 
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(9) 

This equation is diagonally dominant and there- 
fore may be solved through the use of a tridiagonal 
method. The coefficients /?o, 0u an( l A 2 were de- 
termined by numerical experimentation to be 0.16, 
1.00, and 0.50, respectively (ref. 17). These values 
were used initially for all design cases. These coef- 
ficients are discussed in more detail subsequently in 
this paper. 

The design procedure is initiated after an initial 
solution for the starting geometry has been obtained. 
Target pressures are input to begin the design proce- 
dure. All modifications are applied to wing sections 
at the input span stations rather than at the compu- 
tational span stations. If a target pressure distribu- 
tion is desired at a span station that does not corre- 
spond to an input station, a new user-defined wing 
section is generated at that span location through 
interpolation. 

As stated previously, modifications to the wing 
(or airfoil) sections consist of increments to the ge- 
ometry normal to the wing plane. These incremen- 
tal changes are applied to user-defined wing sections 
at each computational wing fine-grid point in the x- 
direction; that is, incremental changes are applied at 
every 1-percent-chord point downstream of the lead- 
ing edge. The leading-edge and trailing-edge points 
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of the wing sections are held fixed since their sur- 
face velocities cannot be accurately computed by this 
analysis code. Boundary conditions are then recalcu- 
lated for the modified wing geometry. Several analy- 
sis iterations, usually 10, are performed on the mod- 
ified geometry and the pressures are again compared 
with the target pressures. This cycle of comparing 
pressure coefficients and modifying the geometry is 
referred to as a design step. The number of design 
steps needed to reach the target pressure distribu- 
tion varies according to the severity of the changes 
required to the initial wing (or airfoil) geometry. This 
is not to say that the target pressure distribution is 
reached for every design case. There are many situ- 
ations in which the design cannot be obtained, such 
as specification of an unrealistic target pressure dis- 
tribution or too severe a change between the target 
and initial pressure distributions. At this time no 
test is used to determine if the designed pressure dis- 
tribution is reasonably close to the target pressure 
distribution. Instead, the user specifies the number 
of design steps to be taken. 

The design code is set up to redesign the up- 
per surface and the lower surface independently, so 
single-surface design (holding the other surface fixed) 
is an option. Since all modifications are made to the 
actual geometry, the geometric modifications can eas- 
ily be constrained by the user if desired. 

Design Test Cases 

Several design test cases were considered to ver- 
ify the compatibility of the design algorithm and the 
analysis code and to demonstrate the capabilities of 
the design algorithm. Two-dimensional cases were 
initially considered, the first being a subcritical tar- 
get pressure distribution and the second a test case 
where there is a significant region of supercritical flow 
and a well-defined shock. A three-dimensional test 
case then followed for a swept and tapered planform. 
Also, two previous design efforts with a cut-and-try 
approach were repeated utilizing the automated com- 
putational design code. For the test cases presented, 
aerodynamic calculations were made on the final de- 
sign geometry and the resulting pressure distribution 
was compared with the target pressure distribution. 

The first two design test cases were two- 
dimensional and used an NACA 0006 airfoil as the 
starting geometry. The target pressure distribution 
for both cases was that of an NACA 0012 airfoil at 
the particular flow conditions. The first case was ex- 
ecuted at a free-stream Mach number of Moo — 0.70 
and a = 1.0°. This case required 30 design steps to 
approximate the target pressure distribution. (See 
fig. 5(a).) Comparisons of the design airfoil ordinates 


with the NACA 0012 ordinates yielded root-mean- 
square (rms) errors of 0.00262 for the upper surface 
and 0.00215 for the lower surface. The rms error was 
much larger than might be expected considering the 
good agreement of the design and target pressures 
shown in figure 5(a). As shown in figure 5(b), the dif- 
ferences in the surface ordinates appear to have been 
the result of an angle-of-attack change. However, 
the leading-edge point (x/c = 0) and the trailing- 
edge point ( x/c = 1.0) were held fixed during the 
design process and angle of attack was not allowed 
to change. It was then conjectured that the cam- 
ber line of the design airfoil differed from that of the 
NACA 0012, which could simulate the target airfoil 
at a lower angle of attack. In figure 6 the camber 
lines of the design and target airfoils are shown. As 
is shown in the figure, the camber line of the design 
airfoil was below that of the NACA 0012. The de- 
sign airfoil was then rotated about the trailing edge 
0.15°, based on the angle of the camber line of the 
design airfoil, and compared with the NACA 0012 
airfoil (fig. 7). The rms error was reduced to 0.00134 
on the upper surface and 0.00123 on the lower sur- 
face. Aerodynamic calculations were made for the 
NACA 0012 at Moo = 0.70 and a = 1.0° and 0.85°; 
these calculations showed no discernible changes in 
the pressure distribution. Based on the informa- 
tion obtained through rotation of the design airfoil 
and by assessment of the effects of small angle-of- 
attack changes on the pressure distribution of the 
NACA 0012, the design geometry obtained for this 
design test case was reasonable. This design case il- 
lustrated a certain lack of sensitivity of the target 
pressure distribution to small changes in the design 
geometry. 

The next case, at Moo = 0.80 and a = 0°, was a 
more difficult design problem in that there is a sig- 
nificant difference in the extent of supersonic flow 
over the NACA 0012 and the NACA 0006 airfoil 
at Mqo = 0.80. Fifty design steps were needed to 
match the target pressure distribution (fig. 8(a)). A 
comparison of the design and NACA 0012 airfoils is 
shown in figure 8(b). The rms error between the 
design and the NACA 0012 airfoil was 0.00163 for 
the upper and lower surfaces. Oscillations which oc- 
curred in the pressure distribution at the base of the 
shock were eliminated after an additional 140 aero- 
dynamic calculations on the final design geometry 
(fig. 9). 

Upon completion of these two-dimensional design 
test cases, a three-dimensional design test case was 
performed. For this case, a swept, tapered planform 
was used for the starting geometry, with the entire 
wing made up of NACA 0006 airfoils (fig. 10). The 
planform had 20° of leading-edge sweep, an aspect 
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ratio of 7.2, and a taper ratio of 0.26. Only airfoil 
sections at two locations were used to define the wing, 
one at the root and one at the tip. The design task 
was to modify the wing section at 77 = 0.665 so that 
the pressure distribution for the NACA 0012 airfoil 
at the design conditions would be obtained at this 
wing station. 

This was a difficult design test case because at 
Moo = 0.90 and a = 0° a strong shock occurred in 
the target pressure distribution at 77 = 0.665, where 
the NACA 0012 airfoil was located (fig. 11). Fig- 
ure 11 shows the design and target pressure coeffi- 
cients after 40 design steps for rj = 0.665 and for 
several other span stations. Oscillations about the 
shock were caused by irregularities in the surface 
slopes and curvature in that region. The pressure 
distributions at span stations other than rj = 0.665, 
although labeled as target pressures, were actually 
known pressures for the desired configuration. As de- 
scribed above, the desired configuration had the same 
planform as the initial geometry with NACA 0006 
sections defining the root and the tip and with an 
NACA 0012 airfoil at r] — 0.665. Only the air- 
foil at 77 = 0.665 was changed with the design 
algorithm. Computational airfoils inboard and out- 
board of 77 = 0.665 were determined through span- 
wise interpolation, which accounts for the oscillations 
in the pressure distributions at span stations other 
than the design station. A comparison of the design 
airfoil ordinates with the NACA 0012 airfoil ordi- 
nates (fig. 12) resulted in an rms error of 0.001191 
for the upper and lower surfaces. The design airfoil 
ordinates were smoothed with a least-squares fit to 
a cubic polynomial with a square root term added 
in the leading-edge region. The smoothed ordinates 
were substituted for the initial NACA 0006 ordinates 
at rj = 0.665 and an analysis run was subsequently 
made on this new wing geometry. Results from the 
analysis run compared well with the desired target 
pressures and are illustrated in figure 13. 

Also considered were two additional test cases 
based on actual designs previously completed by 
means other than automated computational design. 
One was originally addressed by Campbell, Wag- 
goner, and Phillips (ref. 4), who relied on computa- 
tional analyses during the design of a natural laminar 
flow (NLF) airfoil at transonic conditions. The de- 
sign approach was based on the cut-and-try method 
with an airfoil analysis code. The starting point for 
this design was a modified 14-percent-thick NLF air- 
foil developed by Viken (ref. 24). Viken’s original 
airfoil was modified to lower the design lift coefficient 
by reducing the trailing-edge camber. Analysis of the 
initial geometry at M 0 0 = 0.70 and a = —0.193° 
revealed two undesirable features in the pressure 


distribution of this airfoil (fig. 14) that were elimi- 
nated using the cut-and-try design approach. The 
pressure peak at the leading edge of the lower sur- 
face was undesirable since it could cause transition of 
the boundary layer. The second undesirable feature 
was the steep pressure recovery region on the upper 
surface, which could cause the flow to separate. The 
expansion ahead of the shock was not present in the 
results from the full-potential methods used in refer- 
ence 3. The final-design airfoil of Campbell et al. is 
designated the HSNLF-0213. Figure 15 illustrates the 
initial pressure distribution, the target pressure dis- 
tribution (of the final airfoil of ref. 4), and the design 
pressure distribution after 40 design steps using the 
current design code. Although the target and de- 
sign pressure coefficients were reasonably matched, 
there is a significant difference between the design 
airfoil and the HSNLF-0213 airfoil. One factor that 
may have contributed to the discrepancies in this de- 
sign case is a possible path dependency of the design 
algorithm. This was the only design test case at- 
tempted where the initial geometry was thicker than 
the desired final geometry. In order to investigate 
this possible problem, another design case based on 
the HSNLF-0213 airfoil but with a flat plate as the 
initial geometry was then attempted. The results 
after 30 design steps (fig. 16) showed that the de- 
sign pressures again matched the target pressures 
very well, but the agreement between the design air- 
foil and the HSNLF-0213 airfoil was much improved. 
These results indicate that the direction of change 
to the airfoil surface might be critical in the design 
process. 

The second design case attempting to duplicate 
a previous design task was based on modifications 
made by hand to the root of a supercritical wing 
as part of a research study on supercritical wings 
(ref. 25). The wing planform selected for modifica- 
tion had 30° of sweep at the quarter-chord of the 
outer portion of the wing and an aspect ratio of 9.8. 
A leading-edge extension was added at the root to 
increase structural depth (fig. 17). The coordinates 
for the final wing geometry were available and used 
to generate target pressure distributions for wing sec- 
tions inboard of the break. The initial geometry for 
the inboard sections could not be located, so the 
NACA 0006 airfoil was used as the starting geom- 
etry in the region to be designed. The fuselage was 
modeled by an axisymmetric body with a maximum 
radius of 11.7 in. Figure 18 illustrates the target 
and design pressure distributions at Moo = 0.83 and 
a = 2.19° for several span stations after 40 design 
cycles. Oscillations occurred in the pressure distribu- 
tions at all design span stations. Based on the third 
design test case, the design airfoils were smoothed 
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and used as the starting geometry for a second de- 
sign run. After 30 additional design steps the results 
had improved, but oscillations were still present at 
the more outboard stations. The design airfoils were 
again smoothed and used as the starting geometry for 
a third design run. Figure 19 shows the results af- 
ter 30 additional design steps. Good agreement with 
the target pressure coefficients was obtained at the 
inboard stations, but as the shock strength increased 
outboard so did the oscillations in the pressure coeffi- 
cients at the shock. Comparisons of four of the design 
airfoils with the actual airfoils used to generate the 
target pressures are shown in figure 20. Excellent 
agreement is shown for the lower surface, even in the 
cove region. Differences in the upper surfaces vary 
with each span station, with the largest discrepancy 
occurring at the most outboard station. This was 
the most difficult design test case attempted. Super- 
critical airfoils such as those desired in the inboard 
region are very sensitive to small changes in the flow 
conditions, and this sensitivity, as shown in previous 
design cases, complicates the design task. The algo- 
rithm did, however, modify the wing sections in such 
a way as to result in a wing geometry much closer to 
that desired than to the initial wing geometry. 

Design Code Coefficient Selection 

The design algorithm used by the subject code 
contains three coefficients, /?o, an( i f° r which 
the values can be changed by the user. The de- 
fault values have been determined by numerical ex- 
perimentation to be 0.16, 1.00, and 0.50 for /?o, /?i, 
and 02 , respectively (ref. 17). These values, how- 
ever, were not believed to be the optimum values 
for all design cases. A series of design computations 
were performed to study the effects of the values of 
/? 0 , /?i, and /?2 on the design process and to deter- 
mine how these values should be selected for a par- 
ticular design problem. The values of flo, fti, and /?2 
were chosen such that all either increased or de- 
creased by the same factor over the default values. It 
was also found during this series of design computa- 
tions that the value of the subsonic relaxation factor 
u>, which is used by the analysis code to accelerate 
the flow solution, also has an effect on the design re- 
sults. Preliminary guidelines based on these design 
runs for using the design code were determined to be 
as follows: 

1. For design test cases where both the 
initial and final pressure distributions are 
subcritical, use /?o = 0.16, /?i = 1.00, 

/?2 = 1*00, and u = 1.5. 

2. For design test cases where there is a 
significant difference in the amount of 


supersonic flow between the initial and 
final pressure distributions, use /?o = 0.32, 

P 1 = 2.00, /? 2 = 1-00, and tu = 1.8. 

All design cases presented in this paper utilized the 
above guidelines. 

Upon completion of the initial design test cases, 
additional design cases were run for which the values 
of the design coefficients were allowed to change rel- 
ative to one another. Several of these design compu- 
tations were performed with two of the three design 
coefficients equal to zero. The third coefficient var- 
ied for each design run to determine its individual 
effect on the design. In addition, design computa- 
tions were performed for which the values of the de- 
sign coefficients were determined by the previously 
discussed preliminary guidelines, with the exception 
of one coefficient. For this particular coefficient, a 
set of design computations was made in which the 
value of this coefficient increased and decreased from 
the preliminary guideline value. This procedure was 
followed for all three design coefficients. All design 
computations were made at both supersonic and sub- 
sonic conditions with the NACA 0006 airfoil as the 
initial geometry and the target pressure distribution 
of the NACA 0012 airfoil at the design conditions. 

The results of the above design computations 
showed that the coefficient (i\ was more influential 
than the other two coefficients when the flow field 
was primarily supersonic. Also demonstrated was the 
stronger influence of the /?2 coefficient for subsonic 
flow. The coefficient /3q had virtually no effect on the 
design for both supersonic and subsonic conditions. 
These results are in accordance with those obtained 
by Davis (ref. 22) and Campbell and Smith (ref. 26), 
each of whom developed similar design algorithms. 

Concluding Remarks 

A three-dimensional transonic design code has 
been developed that predicts wing contour modifi- 
cations to achieve a specified pressure distribution. 
The analysis portion of the design code is a transonic- 
small-disturbance code applied to the wing, body, 
pods, pylons, and winglet. This analysis code has 
been used extensively on a wide variety of configu- 
rations, a fact which lends confidence to the results 
obtained when coupled with a design algorithm. An 
additional benefit of using this code is its capabil- 
ity for two-dimensional analysis, which means two- 
dimensional design is an option. The design algo- 
rithm employed relates changes in surface velocities 
to changes in the geometry. 

Several two- and three-dimensional design test 
cases were considered to demonstrate the capabilities 
of the design code. The design cases were generally 
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successful, the exception being a three-dimensional 
supercritical-wing design case. The results of the 
design test cases have led to the following conclusions 
and recommendations concerning the design code: 

1. The design code developed during this study 
is operational, as demonstrated by various two- and 
three-dimensional test cases. 

2. Based on initial studies, two sets of values 
for the design coefficients /?o, /?i , and /% and for 
the subsonic relaxation factor have been determined. 
One set of values is for supersonic target pressure 
distributions and the other is for subsonic target 
pressure distributions. 

3. Based on a series of design runs made after the 
completion of the design test cases, it was found that 
for supersonic flow the design coefficient / 3 \ was more 
influential than /?q or /?2 on the design. For subsonic 
flow, /?2 was the most influential. The coefficient 
/?0 had virtually no effect on the design at either 
supersonic or subsonic conditions. 

4. Although the target pressure coefficients are 
achieved, additional design steps may be necessary 
because of the lack of sensitivity of the final wing (or 
airfoil) shape to small changes in the flow conditions. 

5. Further studies need to be conducted on the 
effects of starting with a wing (or airfoil) that is 
thicker than the final-design wing (or airfoil). 

6. In order to glean the full benefits of the de- 
sign algorithm, a systematic study to determine the 
optimum values for /?q, /?i, and /?2 for various flow 
conditions needs to be conducted. 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 
February 5, 1988 
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grid systems 

Figure 1. Computational fine-grid system. (From ref. 1.) 



Figure 2. Crude grid stretching in x- and ^-directions. (From ref. 1.) 
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Figure 3. Wing fine-grid system. (From ref. 1.) 


Design airfoil 


Initial airfoil 



Figure 4. Contour changes to surface geometry. 





x/c 

(a) Comparison of wing pressure distributions. 



(b) Comparison of airfoil geometries. 

Figure 5. Results for design test case with NACA 0006 airfoil as initial geometry and target pressure distribution 
of NACA 0012 airfoil at Moo = 0.70 and a = 1.0°. 
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— Design geometry 
— Target geometry 








x/c 

(a) Comparison of wing pressure distributions. 



(b) Comparison of airfoil geometries. 

Figure 8. Results for design test case with NACA 0006 airfoil as initial geometry and target pressure distribution 
of NACA 0012 at Moo = 0.80 and o = 0°. 
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Figure 11. Results for design test case of swept, tapered 
and o = 0°. 
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Design results 

Targets 

Initial results 
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with modified airfoil at t) = 0.665 at Moo = 0.90 
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Figure 11. Concluded. 
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Design results, smoothed 

Targets 

Initial results 
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n = .665 

Design results, smoothed 

** + + «■* + Targets 
Initial results 



Figure 13. Results for new wing geometry with smoothed airfoil ordinates at tj = 0.665 at Moo = 0.90 and 
a = 0°. 
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n = ,730 n = .874 

i ? Design results, smoothed 

I Targets 
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Figure 13. Concluded. 
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Design results 

Targets 

Initial results 
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(a) Comparison of airfoil pressure distributions. 


Design geometry 

Target geometry 



(b) Comparison of airfoil geometries. 

Figure 15. Results for design case with modified 14-percent-thick airfoil as starting geometry at = 0.70 
and a = -0.193°. 
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(a) Comparison of airfoil pressure distributions. 



(b) Comparison of airfoil geometries. 

Figure 16. Results for design case with modified 14-percent-thick airfoil as starting geometry at M 
and a = -0.193°. 




Figure 17. Planform definition for supercritical-wing design test case. Sweep = 30°; Aspect ratio = 9.8. 
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